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Abstract 



Nonrelativistically forbidden, single-photon transition rates between low ly- 
ing states of the helium atom are rigorously derived within quantum elec- 
trodynamics theory. Equivalence of velocity and length gauges, including 
relativistic corrections is explicitly demonstrated. Numerical calculations of 
matrix elements are performed with the use of high precision variational wave 
functions and compared to former results. 
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The existence of nonrelativistically forbidden transitions in helium, for example between 
the singlet and triplet states, indicates the presence of relativistic effects. The calculation 
of these effects in atoms or ions is a highly nontrivial task. Depending on the magnitude 
of nuclear charge Z one performs various approximations. Here we study light atoms, so 
the expansion in the small parameter Z a is the most appropriate. Forbidden transitions 
have already been studied for many light atoms and especially for helium (for a review see 
|TJ). Historically, the first but approximate calculations of S-P forbidden transitions were 
performed by Elton in ||. Since the dominant part comes from 2 3 Pi and 2 X P\ mixing, he 
included in the calculation only these states. Drake and Dalgarno in |J were the first to 
include higher excited states, which led to much higher precision. Moreover, Drake later 
Q accounted for corrections to S-state wave functions. Although these calculations were 
correct, there was no proof that they are complete. As an example may serve the 2 3 Si — l 1 ^ 
Ml transition. Feinberg and Sucher || derived an effective operator for this transition and 
showed the cancellation of electron-electron terms. However, the calculations of Drake in 
were performed earlier with the implicit assumption, that these terms are absent. In a 
completely different approach based on relativistic many body perturbation theory Johnson 
et al. |l| and Derevianko et al |7j studied forbidden transition in both velocity and length 
gauge. They pointed out the significance of negative energy states. However, not all results 
were in agreement with the nonrelativistic approach based on the Breit hamiltonian. It is 
the purpose of this work to systematically derive matrix elements for forbidden transitions 
in helium within quantum electrodynamics theory. The equivalence of length and velocity 
gauges for El transitions, including relativistic corrections, is explicitly shown. With the 
use of optimized numerical wave functions, the amplitudes and transition probabilities for 
2 3 P 2 - 1%, 2 3 Px - 1%, 2 1 P 1 - 2 3 S U 2 3 S X - 1%, and 3 3 S 1 - 2 3 S X are calculated with high 
precision and compared to former results. 

The nonrelativistic helium atom interacting with the electromagnetic field is described 
by the Schrodinger-Pauli hamiltonian: 

H _ (p 1 -eA) 2 | (p 2 -eA) 2 a Za Za 
2 m 2m r r\ ri 

The single photon transition amplitude T between two eigenstates <j) and ip, in the electric 
dipole approximation is 

T = (0| (Pl+P2) V ) =i{E+- E+) (0|(n + r 2 y\i,) , (2) 
m 

and the transition probability A is 

/ Li Uj \ 

A = 2a\E 4> - E^\TT*- J (>' - -^J . (3) 

In the effective Hamiltonian approach relativistic corrections enter in two ways, as corrections 
to the wave functions 4> and ip and the correction 5] to the current p/m 

f = W sM + <rf*±£ j^jjy am + wh *±^> . (4) 
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The correction to the wave function is given by the Breit hamiltonian. The part responsible 
for singlet-triplet transition is 



5H 



-a — 2 ~3 x Pi ~ "3 x P 2 + 1 — 2 T x (Pi + P 2 
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(5) 



Corrections to the current are given by several time ordered diagrams, shown in Fig. 1. The 
corresponding expression is calculated as follows. The first diagram is 



Sji = u + (p')au(p) 
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where u(p) is a normalized plane wave solution of the free Dirac equation. For considered 
transitions one may leave spin dependent terms only. In position representation it takes a 
form 
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The photon momentum k is of order m(Za) 2 , while r is of order (mZa) 1 . This means 
that e l k ' r can be expanded in powers of k-f. After adding contributions from both electrons 
the (Za) 2 correction takes the form 



Sjt = - — (k ■f 1 )kxa 1 + - — (k ■ r 2 ) k X a 2 . 
2m 2m 



(8) 



The next diagram involves one electron-positron pair and the corresponding expression is 

Ze 2 1 
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(9) 



where A_ is a projection operator into the negative energy subspace and q is a momentum 
exchange between electron and the nucleus. The [Z a) 2 correction from both electrons 
becomes 
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(10) 



The remaining diagrams involve electron-electron terms. The last two are of higher order, 
so they will not be considered here. The expression for diagram 3 can be obtained from 
Eq.(^|) by the replacements — Z a — > a. In this way one obtains 
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where r denotes here r = r*i 2 = f% — r 2 . The (Z a) 2 correction is 



$33 
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2 m 2 r 3 

The expression for diagram 4 is 



— r x (cti - cr 2 ) . 
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(13) 



The term in the second line equals 2 8 %3 and that in the first line has already appeared in 
Eq. (Eh, so it becomes 
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The (Za) 2 correction is 
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and cancels out with that from diagram 3, Eq. ([12]) . The final expression for the relativistic 
correction to the current of order 0(Z a) 2 is the sum of (§) and fllQf ) 



1 - 1 - 1 Za 
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This Sj could be also derived through the Fouldy-Wouythusen transformation of a l e lk ' r , 
however in this way possible electron-electron terms are omitted, which happens to be correct 
for just this case. Having 8j and SH, the transition amplitude T l in (f|) will be transformed 
to the length gauge with the use of identity 



P1+P2 



rn 



i [H, fi + r 2 ] 



(17) 



and the fact the terms in T % proportional to k % do not contribute to the transition rate, as 
it can be seen from Eq. (|3]). After performing simple algebraic transformations the result is 



T*=i(E+-Ei,)U<l>\(r\+r\) 
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(19) 



4 



The first term in Eq. (|18"D corresponds to electric dipole, and the second one to magnetic 
quadrupole transitions. It is worth noting that for electric dipole transitions, as given in 
length gauge, relativistic corrections enter only through corrections to the hamiltonian SH. 

So far, we have considered only forbidden transitions with spin change between S and P 
states, namely 2 3 P 2 — * ^Sq, 2 3 P\ — > l 1 ^ and 2 1 P\ — > 2 3 Si. However, even more forbidden 
Ml transitions 2 3 5i — > l 1 ^ and 3 3 Si — > 2 3 Si arrive at the order O(Za) 3 , so they are not 
described by the expression in Eq. flip]). No second order type of terms contribute and in 
the calculation of 8]m one takes the next corresponding term in the expansion of e lk ' r in 
Eqs.(|,|,0,0) 
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This result agrees with the former one, obtained by Feinberg and Sucher in 
transition between 2 3 S' 1 and l 1 ^ it could be further simplified to 
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k 2 in the above can be replaced by 
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in this way one obtains for SjM another simple expression 
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The analogous expression for the 3 3 Si — 2 3 S' 1 transition reads 
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For Ml 
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We now consider the spin algebra in the calculation of the transition probability, as given 
by Eqs. (|3]) and (|18f) . One summs up over final states and averages out over initial states. 
The appropriate formulas are: 

(25) 
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(27) 
(28) 

(29) 
(30) 
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where s = 0i/2 + 02/2 and the following normalization is utilized: (P l \Pi) = 5^/3. More- 
over, for this calculations one needs two formulas for spin product 



01 - a 2 Y 1 - 



(01 - 2 1 



(01 - 02)* s 2 (01 - o 2 y 

and the following set of formulas for spin traces 

Tr s i = , 
Tr s* s j = 2 5 ij , 
Tr s i s j s k = i e ijk , 
Tr s i s j s k s l = 5 ij 5 kl + 5 jk 8 il . 
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With the help of the above formulas one obtains for transition probabilities (for simplicity 
we put m = 1) the following expressions 
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where k = \AE\, and h % is defined by Eq. (|J). It is worth noting that l P\ — > 3 Si is not only 
El transition but also M2, which has not yet been recognized in the literature. 

Once transition probabilities are expressed in terms of matrix elements between non- 
relativistic wave functions, they can be calculated numerically with high precision. In the 
numerical calculation we follow an approach developed by Korobov ||. The wave function 
is expressed in terms of exponentials 
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The parameters aj,/3j,7j are chosen randomly between some minimal and maximal values, 
which were found by minimization of energy of a specified state. The maximal dimension 
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of this basis set was 600. Lower values were used for checking convergence. The advantage 
of this basis set is simplicity of matrix elements, which are expressed in terms of integral 
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For some more singular matrix elements an additional integral with respect to corresponding 
parameters has to be performed. The disadvantage of this basis set is the necessity of using 
quadruple precision for N > 100. Moreover, the second order terms require more careful 
tuning of parameters due to the singularity of 5H and large mixing of 2 3 P 1 and 2 1 P 1 states. 
These, which involve odd parity intermediate P-states, are much larger than those which 
involve even parity P-states, by approximately three orders of magnitude. It is due to the 
fact that energies of even parity P states lie beyond the ionization level. Most often, these 
small second order terms were neglected in the former calculations. However, they are not 
neglected here. Our numerical results for forbidden transitions between low lying states are 
presented in Table I. 

In the comparison with former work we start with the Ml transition 2 3 S\ — > l 1 ^. This 
transition was measured by Moos and Woodworth in with the result A = 1.10(33) 10~ 4 s _1 
and Gordon Berry from Notre Dame is currently preparing a more precise measurement. 
The first (correct) theoretical result obtained by Drake in || 1.272 • 10~ 4 s _1 is in agreement 
with the experimental value. However, as pointed out by Feinberg and Sucher in ||, Drake 
has not considered electron-electron terms, which happened to cancel out for this transition. 
Later, Johnson et al [|TJ used RMBPT to calculate forbidden transitions for any helium-like 
ions and obtained a result for Z = 2, which is 1.266 ■ 10 _4 s _1 . It differs slightly from 
the result obtained here 1.272426 ■ 10 _4 s _1 , due to inclusion in |1] of some higher order 
terms, while electron correlations were not well accounted for. Moreover there are unknown 
radiative corrections and exchange type of diagrams of order a/(2 7r), the last two in Fig. 
(1), to any of these transitions. Therefore only first 3 digits are physically significant. 
Numerical results are presented with higher precision for the purpose of comparison with 
former results. Next the Ml transition 3 3 Si — > 2 3 S± rate was obtained only by Derevianko 
et al in H. Their result 1.1710 -8 s -1 , disagrees with ours, 6.484690 10~ 9 s -1 . The reason of 
this discrepancy is left unexplained. It may indicate the loss of accuracy of RMBPT due to 
strong numerical cancellation. This discrepancy does not have experimental impact since 
this rate is too small for Z = 2 to be measured. However, calculations should be verified for 
higher Z, where this transition rate grows with Z 10 and becomes measurable at some value 
of Z. The next considered transition is M2: 2 3 P 2 — > l 1 ^. It was first obtained by Drake ||: 
A = 0.327 s -1 , and later by Johnson et al A = 0.3271 s _1 , in agreement with our result 
A = 0.3270326 s -1 . The calculation of the intercombination El transition 2 3 Pi — > l 1 ^ 
was little more elaborate, since it involves infinite summation over intermediate states. In 
former works the second term in Eq. ( p7[ ) involving even parity P-states was neglected. 
Indeed, calculations show it is smaller than 1%. The first complete result by Drake [|J is 
A = 176.4 s -1 . RMBPT calculations of Johnson et al |T| including negative energy states 
is A = 175.7 s _1 and our result A = 177.5771 s _1 agrees within 1%. The last transition 
2 : Pi — ► 2 3 S\ is a sum of El and M2. The result A = 1.55 s _1 obtained by Drake includes 
only El transition. Our result is A = 1.548945 s _1 and the magnetic transition happened to 
be negligible 0.000019 due to small energy splitting. 
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In summary, we have presented a rigorous derivation of rates for nonrelativistically for- 
bidden transitions. We demonstrated equivalence of length and velocity gauges including 
relativistic correction for forbidden transitions. We confirmed the commonly used fact that 
in the length gauge relativistic corrections enters only through corrections to wave function 
as given by Breit hamiltonian. We verified that M2 2 1 P\ — > 2 3 S*i transition is much smaller 
than El, which was implicitly assumed in former works. Our numerical calculations using 
simple exponential functions confirmed former results with the exception of 3 3 Si — > 2 3 S± 
transition, where our result is approximately twice smaller than of 0. 
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FIGURES 




FIG. 1. Time ordered diagrams for corrections to the current. Dashed line is a Coulomb photon, 
the wavy line is the transverse photon. 
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TABLES 



TABLE I. Transition rates in helium in units s 1 , [n] = 10 n 



transition 



AE in atomic units 



rate A 



E1+M2: 
El: 
M2: 
Ml: 
Ml: 



2 1 P± 
2 3 Pi 

2 3 P 2 
2 3 S 1 
3 3 5i 



2 3 ^ 
1% 
1% 
1% 

2 3 ^ 



0.0513862917 
0.7705606863 
0.7705606863 
0.7284949988 
0.1065403108 



1.548945 

1.775771[2] 

3.270326[-l] 

1.272426[-4] 

6.484690[-9] 
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